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DMPK Equation for Transmission Eigenvalues in Metallic Carbon Nanotubes 
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The Dorokhov-Mello-Pereyra-Kumar (DMPK) equation for transmission eigenvalues is de- 
rived for metallic carbon nanotubes with several conducting channels when the potential range 
of scatterers is larger than the lattice constant. With increasing system length L, the system 
approaches a fixed point, where only one channel is perfectly conducting and other channels are 
completely closed. The asymptotic behavior of the conductance in the long-L regime is investi- 
gated on the basis of the DMPK equation. It is shown that the length scale for the exponential 
decay of the typical conductance is reduced due to the presence of the perfectly conducting 
channel. If a magnetic field is applied, the system falls into the unitary class. It is pointed out 
that this transition is characterized by the disappearance of the perfectly conducting channel 
and the increase in decay length for the typical conductance. 



Carbon nanotubes (CNs) are micrometer-long hollow 
cylinders with nanometer-scale radii. 1 We can regard 
CNs as a new class of quasi-one-dimensional quantum 
wires. Recent experiments have enabled the electrical 
transport measurement of individual single- wall CNs. 2-4 
We expect that the unique topological structures of CNs 
induce unusual electron transport properties, which have 
not been observed in ordinary quantum wire systems. 
Ando and Nakanishi 5,6 studied electron scattering in 
CNs in the case where only one conducting channel is 
present and showed the absence of backward scattering 
when the potential range of scatterers is larger than the 
lattice constant. 

Recently, Ando and Suzuura 7 studied metallic CNs 
with long-range impurity potential in the case where 
the Fermi level lies in several conducting channels and 
showed within a k ■ p scheme that one perfectly con- 
ducting channel is present. They also performed numer- 
ical simulations based on a tight-binding model. Their 
numerical result suggests that, with increasing system 
length L, the system approaches a fixed point, where only 
one channel is perfectly conducting and other channels 
are completely closed. This denotes that the conductance 
G decreases towards the quantized value Go governed by 
the perfectly conducting channel. Due to the spin and 
valley multiplicity, the quantized conductance is given 
by Go = 4e 2 //i. Takane and Wakabayashi 8 studied the 
same problem by adopting a random-matrix theory and 
showed that the fixed point is stabilized by the antilo- 
calization effect. They also showed that G in the pres- 
ence of electron decoherence behaves as G oc L^/L when 
' ^ L, where I and are the mean free path and 

the phase coherence length, respectively. These results 
indicate that the transport properties of CNs with sev- 
eral conducting channels are very different from those 
observed in ordinary quantum wire systems. However, 
our understanding of this problem is not sufficient at the 
present stage. Indeed, we do not know how the system 
approaches the fixed point with increasing system length. 

To answer such a question, we must investigate the 



scattering problem in CNs from a statistical viewpoint. 
Generally, electron scattering in a quantum wire system 
with TV channels is described by the scattering matrix s, 
which is a 2N x 27V matrix given by 9 



f 



(1) 



with N x N reflection matrices r and r' (reflection from 
left to left and from right to right) and transmission ma- 
trices t and t' (transmission from left to right and from 
right to left). The Hermitian matrices t't'^ and ttf have 
the same set of eigenvalues T\,Ti,--- ,Tjy. The trans- 
port properties of a sample are determined by {T a }. For 
example, the conductance is given by G = (2e 2 /h)T 
with T = J2 a=1 T a , where the factor two corresponds to 
the electron spin. In some cases, it is convenient to use 
A a = (1 — T a )/T a instead of T a . Let us call T a and A a the 
transmission eigenvalues. If the probability distribution 
of the transmission eigenvalues is obtained as a function 
of L, we can completely describe the statistical proper- 
ties of G. In ordinary quantum wire systems, the evolu- 
tion of the distribution function with increasing L is de- 
scribed by the Dorokhov-Mcllo-Pereyra-Kumar (DMPK) 
equation. 9-11 We can obtain useful information from it 
in considering the statistical properties of G. 

In this letter, we derive the DMPK equation for the 
transmission eigenvalues in metallic CNs with long-range 
impurity potential and, on the basis of the resulting 
DMPK equation, investigate the asymptotic behavior of 
the conductance in the regime where the system length 
L is much longer than the mean free path I. In deriv- 
ing the DMPK equation, we adopt the random-matrix 
model presented by Takane and Wakabayashi 8 and con- 
sider the unique scattering symmetry observed in CNs. 
With increasing L, the conductance decreases towards 
the quantized value Go- We calculate several statistical 
averages which characterize the asymptotic behavior of 
the conductance. It is shown that the length scale for the 
exponential decay of the typical conductance is reduced 
due to the presence of the perfectly conducting channel. 
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If a magnetic field is applied, the system falls into the 
unitary class. We point out that this transition is charac- 
terized by the disappearance of the perfectly conducting 
channel and the increase in decay length. 

The number of conducting channels, N, in a metallic 
CN is odd, so we set N = 2m + 1 (m: integer). Ando 
and Suzuura 7 studied the scattering problem when the 
potential range of scatterers is larger than the lattice 
constant and obtained the symmetry relations for the 
elements of s. Their result reads 



H = t', 

r = —r and 



(2) 
(3) 



Equation (2) holds for arbitrary systems with the time- 
reversal symmetry, while eq. (3) characterizes the pecu- 
liarity of our system. Equation (3) directly results in the 
absence of backscattering, r aa — 0. 5,6 Due to the valley 
multiplicity, the conductance is given by G = 4(e 2 //i)I\ 
We introduce the transfer matrix M, which can be ex- 
pressed in terms of the elements of s, 



M = 



r't' 



(4) 



We employ the parameterization' 



8, 12 



M 








77*77) 



(5) 



[l + rfq')* 77 
77* (14 

with 6 — ih, where h is an arbitrary N x N Hermitian 
matrix and 77 is an arbitrary NxN complex matrix. They 
satisfy #* = —9 and t r/ = —77. Equation (5) ensures both 
the flux conservation and the unique scattering symme- 
try. In terms of 9 and 77, we can express the elements of 
s as 

f = c e (l+?7?7 t )-3 7 
t' = (l + rfr])-h- \ 
r = —(1 + rfrj)~^rf, 
r' = c 6 t](1 +77 t 77)"3e"' 



(6) 
(7) 
(8) 
(9) 



We introduce the transfer matrix Mi for a system 
of length L, which is parameterized by 9\ and r\\ with 
cq. (5). Let t'i be the corresponding transmission matrix, 
and T a and A a be the eigenvalues of t'\t'\ and 77*771, re- 
spectively. We see from eq. (7) that T a = l/(A a + l). Due 
to the skew-symmetric nature of 771, one of A a is equal 
to and all the other eigenvalues are twofold degener- 
ate. This indicates that one of T a is unity, resulting in 
the presence of one perfectly conducting channel, 7 and 
all the other eigenvalues are twofold degenerate. This is 
a peculiar nature of CNs with long-range impurity po- 
tential. We introduce a unitary matrix v, reducing t'\t'\ 
to the diagonal form by the transformation 



,T N ). 



(10) 



vt\t'W = diag(Ti,T 2 , 
We find that 

vt\t'y = (l + ^JTi) -1 , (11) 

where 771 = v*r]\v^. Note that '771 = —771. Without loss of 
generality, we assume that the Nth. channel is perfectly 



conducting (i.e., Ajy = 0) and express the NxN skew- 
symmetric matrix 771 as 

£ \ 

7: = 1 _£ ; > (12) 

•■• / 

where £ is an m x m diagonal matrix given by 

£ = diag V^2; ' ' ' . V^) • (13) 

We easily see that 

77}?7i = diag(Ai, A 2 , • • • , A m , Ai, A 2 , • • • ,A m ,0). (14) 

The above argument leads to T a = T a+m = 1/(1 + A a ) 
for 1 < a < m and Tjy = 1. We shall use the notation, 
a = a + m for 1 < a < m and a = a — m for m + 1 < 
a < 2m. Using a, we can rewrite the degeneracy relation 
T a = T a+m (1 < a < m) as T a = T a (1 < a < N - 1). 

To study the evolution of T a with increasing L, wc 
attach a small segment of length Sx to the right-hand 
side of the system. If the transfer matrix for the small 
segment is SM, the transfer matrix M 2 for the combined 
system is given by 

M 2 = 6MM 1 . (15) 

Let t'i be the corresponding transmission matrix. The 
eigenvalues T' a of £' 2 t' 2 are slightly different from T a , so 
we express T' a = T a + ST a . To obtain the DMPK equa- 
tion for the distribution function P(Ti,T 2 ,--- ,T m ,L), 
we must calculate (ST a ) ss and (ST a 5Tb) ss , where (• • • ) ss 
denotes the ensemble average with respect to the small 
segment. The DMPK equation is given by 9 

I -fa-'t^(-««»-P 

1 m B \ 

+ 2^W b {STaSTb)ssP J- (16) 

We apply the parameterization given in eq. (5) to 6M. 
Following Mello and Tomsovic, 12 we assume that (9) ss = 

(r]) ss = 0, (9 ab f] cd ) ss = (6abVcd)sa = (VabVcd)ss = 0, (9 2 + 

W*) ss = and 

(9ab9 c d)ss = K abcd , (17) 

{9ab9* d ) ss = K a b,cd> (18) 

(Vabvld)™ = K lb,cd- (19) 

We shall take the weak-scattering limit, 12 where a mo- 
ment higher than the second plays no role, so that SM 
can be expressed in the simple form 



5M 



l + i 
77* 



77 

1 + 9* 



(20) 



We can relate t' 2 and t'i by eq. (15). From eqs. (4), (15) 
and (20), we obtain 

t'i = t\ (1 - 9* - 77V'! + (0*) 2 ) • (21) 

We calculate (ST a ) ss and ((5T a <5T;,) ss by perturbation 
theory. To second order in perturbation theory, we obtain 

w ab w ba 



ST a = w aa + 



6(#o) 



(22) 
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where w a b is an element of the Hermitian matrix t' it'\ ~ 
t'it'\ in the basis where t' \t'\ is diagonal. The Hermitian 
matrix is given by 

w = vt\ (V) 2 + rfr'Wlr) - r'\r] - ijVi) t'\vl (23) 

We here take the weak-scattering limit; 12 Sx approaches 
and at the same time, we let scattering in the small 
segment become infinitely weak, so that 
-l . 



lim (Sx) < 6 d = a l ab d 

ox — '0 



(24) 



are finite quantities, while (Sx) 1 times a moment higher 
than the second vanishes. Accordingly, we calculate 



f a = lim (Sx) 1 (ST a ) ss , 
f ab = lim (Sx)- 1 (ST a ST b ) s 



k—>0 



(25) 
(26) 



instead of (ST a ) ss and (ST a ST b ) ss . In accordance with 
9^ = —9 and *r/ = —77, we adopt the model 8 in which 

<b,cd is s iven b y 



^ab.cd = lim 1 (0 ab 9 cd ) ss = -S ad 8 bc a' ab , 



• 



a ab,cd 



lim .(Sx) 1 {9 ab 9* d ) ss = S ac 5 bd a' ab , 



5x— *0 



(27) 
(28) 



a lb,cd = }^ n ( Sx ) 1 (Vabvld)ss = ( S adSbc ~ S ac S bd ) <T ab . 

(29) 



»0 



Here, a ab represents the average reflection coefficient per 
unit length from 6 to a, while a' ab represents the average 
transmission coefficient. They satisfy Eb ^ — Eb^ab 
since (9 2 + ?7^)ss = 0. To proceed, we must assume that 
o ab is expressed in a simple form. We adopt the simplest 
choice 8 

aab = w=Tji' (30) 

where I is the mean free path. This ensures the absence 
of backscattering (i.e., a aa = 0). 5,6 Equation (30) leads 

to Eb^ab = Eb<b = 

In calculating f a , we substitute the matrix elements 
of w into eq. (22) and then take the weak-scattering 
limit after the ensemble average with respect to the small 
segment. The ensemble average is carried out based on 
eqs.(27)-(30). Employing eqs. (7), (9) and (11)-(13), we 
obtain /at = and 



f =_Zk 

J a 7 



T, 



(N 1)1 



A" 



T a + T b - 2T a T b 



b=l 
(b#a,a) 



T a -T b 



T 

- 1 a 



+ 



I (N-l)l 



T a + 2 £ 



T a + T b - 2T a T b 



(31) 



b=l_ 
(b^a,a) 



for 1 < a < N — 1. Similarly, we obtain f a N = for 
1 < a < N and 

2 



fab = (Sa,b + Sg,,b) 



jT^(l - T a ) 



(32) 



(N-l)V 

for 1 < a, b < N - 1. Substituting eqs. (31) and (32) into 
eq. (16), we obtain the evolution equation 



dP(T u - 



d 



ds 



a=l 



-faP(T u ---,T m ,s) 



2dT a faaP{Tl '" 



(33) 



where s = L/lis the normalized system length. Although 
we have treated T a as independent variables, it is more 
convenient to use A a = (1 — T a )/T a in some cases. Upon a 
change of variables from T a to A a , we obtain the DMPK 
equation for P(\\, ■ ■ ■ ,X m ,s), 



dP(X 1 , ■■■ ,X m ,s) 



1 



ds 



d\ a 



N 

\ a (l-\a)J 



d 

d\ a V 



^ P(X 1 ,--- ,X m ,s) 



where 



j =ri A c x n n \^- x »\ 

c=l b=l a=b+l 



(34) 



(35) 



Equations (33) and (34) are the central results of the 
present letter. 

From eq. (33), we can derive the evolution equation 
for an arbitrary function F(T\,--- ,T m ). 13 Multiplying 
both sides of eq. (33) by F and integrating over {T a } 7 we 
obtain 



9(F) 
ds 



a=l \ 

Tg 

N-l 



+E 

a=l 



T a 



-Ta + 2 J2 



T a +T b - 2T a T b 



b=l 

(b#a) 



N-l 



T 2 a (l~T a ) 



T a -n 

d 2 F\ 



OF 

df a 



dT? 



(36) 



where the ensemble average (• • • ) is defined by 

/m 
Y[4T a ---P(T 1 ,---,T m ,8). (37) 
a=l 

As an example, we derive the evolution equations for T p , 
where T = 1 + 2J27=i T a- Substituting F = TP into 
eq. (36), we obtain 

d(TP) 



(N-l)- 



ds 



= -p(TP +1 )+p(T 2 TP- 1 ) 

+ 2p(p-l){(T 2 -T 3 )TP- 2 ), (38) 

where T q = 1 + 2 E^Li Equation (38) is in agreement 
with the previously reported result. 8 

Now we consider the asymptotic behavior of the con- 
ductance G in the long-i regime based on the DMPK 
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equation. - Wc introduce a new set of variables x a , re- 
lated to A a by 



A a = sinh x a , 



(39) 



where x a > 0. The perfectly conducting channel is 
characterized by xjy = 0. If we make a change of 
variables from A a to x a , the evolution equation for 
P(xi, X2, ■ ■ ■ ,x m ,s) becomes 



dP 



4(iV- 1) 



E 



d 
dx a 



dP 

dx a 



P dx^j 



(40) 



with 



Q, = -ln(j JJsinli2a; Y (41) 

^ a=l ' 

In the limit of s/N = L/Nl — > oo, we expect that 
the variables x a (1 < a < m) are much larger than 
unity and are widely separated. We assume that 1 <C 
x\ <C X2 <C • • • <C £ TO . Thus, we can approximate that 
Q w — 2^^ =1 (4a — l)x a + constant. Substituting this 
into eq. (40), we obtain 



dP 



1 



V 



d 



dP 

dx a 



4 (JV - 1) ^ &c„ 

v ' a— 1 

The solution of the above equation is 



2(4a-l)P). (42) 



P(x\ , • • • , x m , s) 
where 7 = 2 (AT — 1) and 



— YU 



0=1 



7 ; 



4a- 1 



(43) 



(44) 



We observe that each eigenvalue obeys a Gaussian. 
From the resulting distribution, we obtain (x a+ \ — x a ) = 
4s/7 and ^/Var[a; a ] = \fs~py, where VarL4] = (A 2 ) — 
(A) 2 . This denotes that the variables are widely sepa- 
rated as long as s/7 = L/jl ^> 1. Note that (xi) = 3s/7. 
Thus, when 3s/7 = 3T/7Z ^> 1, the deviation of the con- 
ductance from the quantized value Go = 4e 2 //i is domi- 
nated by x x as SG = G - G ~ G x 8c _2:El . Taking the 
ensemble average, we obtain 



-(hi(SG/G )) = 



Var pn(JG/Go)] 



N — 1 
2 

TV - 1 



(45) 



(46) 



We estimate the decay length £ of the typical conduc- 
tance by identifying exp(ln(<5G/G )} = cxp(— 2£/£). We 
find that 



1)1 = -ml. 

' 3 



(47) 



The decay length corresponds to the localization length 
in ordinary quantum wire systems. However, wave func- 
tions in our system are not localized due to the presence 
of the one perfectly conducting channel, so we do not 



call £ the localization length. The average of SG/Go is 
calculated as 

(<5G/G ) = [°° dxe- 2x e-*( x - 3s ^ 2 

V 2?rs Jo 



2 L_ 

= 8e ; . 



(48) 

This result denotes that (SG/Go) decays more slowly 
than cxp(ln((5G/Go)). Equation (48) is also obtained 
from eq. (38). If wc set p = 1, eq. (38) is reduced to 

(^i)^ = -<r 2 ) + <r 2 >. (49) 

Since the deviation of the conductance is dominated by 
x\ when L/jl 3> 1, we can approximate that T ~ l + 2Ti 
and T 2 ~ 1+2T 2 . Substituting them into the above equa- 
tion and retaining the lowest order terms with respect to 
Ti, we obtain 2{N -l)d(T t ) /ds = -4(Ti). This is easily 
solved as (Ti) oc c 1 , which is equivalent to eq. (48). 

We compare our result with that for the case where the 
number of conducting channels is even (i.e., N = 2m). It 
should be emphasized that there is no perfectly conduct- 
ing channel in the even-channel case. 7 Although the even- 
channel case is not realized in CNs, the comparison pro- 
vides us physical insight into the role of the perfectly con- 
ducting channel. If we adapt our argument to the even- 
channel case, it is reduced to that for the ordinary sym- 
plectic class 9 with m channels. This is consistent with 
Suzuura and Ando's argument 14 that a two-dimensional 
graphite sheet belongs to the symplectic class when the 
potential range of scatterers is larger than the lattice 
constant. In the symplectic class with m channels, the 
localization length is £3=4 = (4m — 2)1. We observe that 
the decay length £ is a factor 1/3 smaller than £,p = i in 
the large-m limit. This reduction should be attributed 
to the presence of the perfectly conducting channel, for 
which a; at = 0. The repulsion between xn — and the 
dominant eigenvalue x\ tends to increase x\, resulting in 
the reduction of £. 

Thus far, we have treated CNs in the absence of a mag- 
netic field. If a magnetic field is applied, the symmetry 
relations given in eqs. (2) and (3) no longer hold. Con- 
sequently, the system falls into the unitary class. This 
transition is characterized by the disappearance of the 
perfectly conducting channel and the increase in decay 
length (£ -> ^ =2 = 2NI). 

In summary, we have studied the quantum electron 
transport in metallic carbon nanotubes with several con- 
ducting channels when the potential range of scatter- 
ers is larger than the lattice constant. Considering the 
unique scattering symmetry observed in carbon nan- 
otubes, we derived the DMPK equation for the distri- 
bution of the transmission eigenvalues. With increasing 
system length L, the conductance decreases towards the 
quantized value Go = 4e 2 jh. We calculated several sta- 
tistical averages which characterize the asymptotic be- 
havior of the conductance, based on the DMPK equa- 
tion. It is shown that the length scale for the exponential 
decay of the typical conductance is reduced due to the 
presence of the perfectly conducting channel. The role of 
a magnetic field is discussed briefly. 
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